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ABSTRACT 

Numerical hydrodynamical simulations have proven a successful means of reproduc- 
ing many of the statistical properties of the Lya forest as measured in high redshift 
quasar spectra. The source of ionization of the Intergalactic Medium (IGM), however, 
remains unknown. We investigate how the Lya forest may be used to probe the nature 
of the sources. We show that the attenuation of Lyman continuum photons by the IGM 
depends sensitively on the emissivity of the sources, permitting a strong constraint 
to be set on the required emissivity to match the measured values of the mean IGM 
Lya optical depth. We find that, within the observational errors, QSO sources alone 
are able to account for the required UV background at z ^ 4. By contrast, the emis- 
sivity of Lyman Break Galaxies must decline sharply with redshift, compared with 
the estimated emissivity at z « 3, so as not to over-produce the UV background and 
drive the mean Lya optical depth to too low values. We also investigate the effect of 
fluctuations in the UV background, as would arise if QSOs dominated. To this end, we 
derive the distribution function of the background radiation field produced by discrete 
sources in an infinite universe, including the effects of attenuation by an intervening 
absorbing medium. We show that for z £ 5, the fluctuations significantly boost the 
mean Lya optical depth, and so increase the estimate for the mean ionization rate 
required to match the measured mean Lya optical depths. The fluctuations will also 
result in large spatial correlations in the ionization level of the IGM. We show that 
the large mean Lya optical depth measured at z ~ 6 suggests such large correlations 
will be present if QSOs dominate the UV background. A secondary, smaller effect of 
the UV background fluctuations is a distortion of the pixel flux distribution. While 
the effect on the distribution may be too small to detect with existing telescopes, it 
may be measurable with the extremely large telescopes planned for the future. We also 
show that if QSOs dominate the UV background at z w 6, then they will be sufficient 
in number to rejuvenate the ionization of a previously ionized IGM if it has not yet 
fully recombined. 
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1 INTRODUCTION 



Numerical simulations of structure formation in the universe 
incorporating hydrodynamics in Cold Dark Matter (CDM) 
dominated cosmologies have proven very successful in repro- 
ducing the statistical properties of the Lya forest as mea- 
sured in high redshift Quasi-Stellar Object (QSO) spectra 
(Cen et al. 1994; Zhang, Anninos & Norman 1995; Hern- 
quist et al. 1996; Zhang et al. 1997; Bond & Wadsley 1997; 
Theuns, Leonard & Efstathiou 1998a). The high level of sta- 
tistical agreement, as high as a few percent in the cumulative 
spectral flux distributions (Meiksin et al. 2001), suggests 



that the models are capturing the essential physical nature 
of these highly ionized absorbing structures. 

Still unknown is the nature of the sources of photoion- 
ization. The most likely candidates are photoionization by 
Quasi-Stellar Objects (QSOs) and stars in young forming 
galaxies or globular clusters. Estimates of the UV back- 
ground required to reproduce the measured Lya optical 
depths of the Lya forest while staying within the nucleosyn- 
thesis bounds of the cosmic baryon density agree within a 
factor of a few with estimates of the QSO contribution to the 
UV background (Meiksin et al. 2001). This suggests that 
QSO sources contribute significantly to the UV background, 
and may even dominate over all other sources. A contribu- 
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tion will also arise from star forming galaxies at z > 3, from 
which Lyman-continuum photons have been detected (Stei- 
del, Pettini & Adelberger 2001). The total contribution from 
QSOs, or Active Galactic Nuclei (AGN) more generally, and 
star-forming galaxies, however, remains uncertain because of 
the flux limitations of the surveys that detect them at these 
high redshifts. Generally only the brightest sources enter the 
surveys, leaving the possibility of a large contribution that 
has gone undetected. Distinguishing the relative contribu- 
tions of young stars vs AGN is thus not possible based soley 
on the average detected contributions of each to the UV 
metagalactic background. 

The discreteness of the sources presents an opportu- 
nity for distinguishing between the possible natures of the 
dominant sources. If star forming galaxies or globular star 
clusters contribute a comparable amount or more of ion- 
izing photons to the metagalactic background as do QSO 
sources, then they must have a much higher spatial den- 
sity than QSOs because of their relatively weaker individ- 
ual luminosities. The resulting UV background would then 
be fairly homogeneous, except for local modulations arising 
from the large scale structure of their spatial distribution. 
By contrast, large fluctuations in the UV background are ex- 
pected at high redshifts (z £ 4) if QSO sources dominate the 
background because of the rapidly growing optical depth of 
the Intergalactic Medium (IGM) to Lyman continuum pho- 
tons and the rapidly diminishing number density of QSO 
sources (Zuo 1992a). The fluctuations in the QSO contri- 
bution will become particularly pronounced at z > 5, where 
the mean Lya optical depth exceeds unity (Fan et al. 2001b; 
Becker et al. 2001), as we show here. 

The role of fluctuations in the UV background on the 
absorption properties of the Lya forest have been largely 
unexplored (but see Croft et al. 2002) . In this paper, we in- 
vestigate the effect of the fluctuations on the 1-pt flux distri- 
bution of the Lya forest. We examine the effect of the back- 
ground fluctuations on the power spectrum of the Lya forest 
in a companion paper (Meiksin & White 2003) . We base the 
contributions of the QSOs to the UV background and its 
fluctuations on the QSO luminosity functions obtained by 
the 2dF (Boyle et al. 2000) and the SDSS (Fan et al. 2001a, 
2001b, 2003). 

Because of the large computational overhead in incor- 
porating proper hydrodynamics into an N-body code, we 
perform simulations of the Lya forest using a pure gravita- 
tional Particle- Mesh N-body code (PM simulations). These 
simulations have been shown to agree well with the full hy- 
drodynamics computations (Meiksin & White 2001), and 
have the advantages of permitting both larger simulations 
to be run as well as greater numbers of them. Ultimately, a 
detailed comparison with the data will require solving the 
full set of hydrodynamics equations. 

This paper is organized as follows: In section 2, we dis- 
cuss our simulations of the Lya forest and model assump- 
tions; estimates of the metagalactic emissivity are given in 
section 3; we derive the distribution function for the UV 
background fluctuations in section 4; in section 5 we discuss 
the pixel flux distribution; we finally summarise our results 
in section 6. 



2 THE SIMULATIONS AND MODELS 
2.1 Particle- mesh simulations 

We generate spectra of the Lya forest from pure DM simu- 
lations, mimicking the temperature of the gas using a poly- 
tropic equation of state. Comparisons with full hydrodynam- 
ics simulations show that simulations recover the principal 
spectral properties of the Lya forest to an accuracy of 10% or 
better in the cumulative distributions of flux, wavelet coeffi- 
cients, H i column density, and Doppler parameter (Meiksin 
& White 2001). 

The evolution of the (dark) matter is computed using a 
parallel implementation of a PM code. The computational 
volume is chosen to be periodic with side length 1, and we 
take as our time coordinate the log of the scale factor, log a, 
where a = (1 + z) _1 . Velocities are measured in units of 
the expansion velocity across the box aHLho* of comoving 
side Lbox, where H is the Hubble parameter fixed by the 
Friedmann equation 

+a > + <$) a, 

in terms of the Hubble constant Ho and the densities in mat- 
ter (Qm), cosmological constant (Qa) and curvature (fix), 
in units of the critical density. The density is defined from 
the particle positions by assigning particles to a regular 
cartesian grid using the Cloud-In-Cell (CIC) scheme (Hock- 
ney & Eastwood (1988)). The Fourier Transform is taken 
and the force computed using the kernel k/k 2 . We use a sec- 
ond order leap frog method to integrate the equations. The 
relevant positions are predicted at a half time step and used 
to calculate the accelerations which modify the velocities. 
The time step is dynamically chosen as a small fraction of 
the cell crossing time with a maximum size of A log a = 2%. 

The initial conditions were created by displacing the 
particles from a uniform grid using the Zel'dovich approx- 
imation. The initial conditions were Gaussian, however to 
avoid large fluctuations between different realizations we 
only accepted those whose low-fc modes closely tracked the 
"mean" . 

Given a set of final particle positions and velocities, 
we compute the spectra as follows. First the density and 
density-weighted line-of-sight velocity are computed on a 
grid (using CIC interpolation as above) and Gaussian 
smoothed using FFT techniques in order to adequately sam- 
ple the velocity field. This forms the fundamental data set. 
A regular grid of 32 x 32 sightlines is drawn through the box, 
parallel to the box sides. Along each sightline we integrate 
(in real space) to find t(u) at a given velocity u. Specifically 
we define 

1 2 



t(u) = / dxA(x) 



p(x) 



T(x) 



-(u-u y/b* 



(2) 



where uq — xaHL\> ox + vi os and b = ^/2fcsT '/ran is the 
Doppler parameter, where mn is the mass of a hydrogen 
atom. The flux at velocity u is exp(— r). The integration 
variable x indicates the distance along the box in terms 
of the expansion velocity across the box. In terms of the 
baryon density parameter fif,, hydrogen baryonic mass frac- 
tion X, the metagalactic photoionization rate F_i2, in units 
of 10~ 12 photons s _1 , and the number f e of electrons per 
hydrogen nucleus, the coefficient A is given by: 
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A(x) = 18.Qf e ^y (^) ri! 2 (x)L box (l + z ) 5 , (3) 

where Lbox is in Mpc and b in eq. (2) is in km s _1 . The ex- 
plicit x dependence of A accounts for any spatial fluctuations 
in the background ionization rate r_i2(x). The electron frac- 
tion is given by / e = 1 + (1 - X)/{2X) « 1.16 (X = 0.76) 
for fully ionized helium, and f e = 1 + (1 - X)/(4X) « 1.08 
if the helium is only singly ionized, a possibility at the high 
redshifts we consider. We shall adopt / e = 1.1 for simplic- 
ity, noting that a different value for / e is effectively absorbed 
into the ionization rate. 

In the absence of UV background fluctuations, A is 
constant and is iteratively adjusted to obtain a predefined 
f a = — log (exp( — t)}. To account for a fluctuating UV back- 
ground field, we replace A by A_i2/[rs,-i2 + f qso,-i2(1 + 
<5qso)], where Fqso,-12 is the average photoionization rate 
produced by the QSO sources, in units of 10 _12 s _1 , A-12 
is the value of A for r_i2 = 1, and <5qso accounts for the 
spatial fluctuations in the QSO contribution. The values for 
<5qso are drawn randomly from the distribution functions 
for the background intensity derived below. We shall con- 
sider two cases for the sources of the UV background. In the 
first, we assume the only sources are QSOs, so that Ts = 0. 
In the second, we shall allow for a smooth component Ts to 
the UV background. This component will partly be present 
due to the re-radiation of ionizing photons through recom- 
bination in the Lya forest itself, at a level by as much as 
(a a — aB)/ctA ~ 40% of the QSO contribution (Meiksin & 
Madau 1993; Haardt & Madau 1996). Because of structure 
in the Lya forest, the re-radiated component will also fluc- 
tuate. We neglect this contribution to the fluctuations since 
they likely will be less than the direct contribution from 
QSOs. Other possible contributions to the ionizing back- 
ground include galaxies, like Lyman Break Galaxies, and 
smaller stellar systems like intergalactic star clusters. In this 
case, the value of Fs is iterated until agreement with f a is 
found. 

In evaluating Eq. (2) we assume a power-law relation 
between the gas density and temperature 

r = T„(|y . (4) 

In practice, To and 7 will depend on the reionization history 
of the IGM, but we shall fix T = 2 x 10 4 K and 7 = 1.5 
throughout for most of our models. The value we use for To 
is larger than used by Meiksin & White (2001), to account 
for the broad lines measured over the full range of line-centre 
optical depths near z — 3.5 (Meiksin et al. 2001). An effort 
to calibrate the evolution of To and 7 was undertaken by 
Schaye et al. (2000). The origin of the extra broadening is 
unknown. A plausible explanation is late He 11 reionization, 
but it may be due to turbulence driven by galactic (or star 
cluster) winds, in which case it may well be present to red- 
shifts as high as 6 or even higher. Since we are not con- 
cerned with detailed matching to the absorption statistics 
of particular spectra in this paper, we have not included a 
complicated evolutionary scenario for To and 7. Instead we 
simply consider a spread in To by also considering cases with 
T = 10 4 K , noting that the equation of state is only an ap- 
proximate description of the IGM temperature and is in any 



case valid only for overdense regions (Meiksin 1994; Theuns 
et al. 1998b). 

2.2 Cosmological model 

We consider a flat cosmological model with Qm = 0.30, 
n A = 0.70, n b = 0.041, h = Jfo/lOOkms- 1 Mpc" 1 = 0.70, 
and slope of the primordial density perturbation power spec- 
trum n = 1.05. The model is COBE normalized, so that the 
fluctuation normalization at 2 = in a sphere of radius 8h~ 
Mpc is cr 8h -i = 0.970, and the fluctuation normalization at 
2 = 3 spherically filtered on the scale of the Jeans length 
for T = 2 x 10 4 K gas is aj(z = 3) = 1.59. The model is 
consistent with current constraints from the CMB (Stom- 
por et al. 2001; Netterfield et al. 2002; Pryke et al. 2002), 
Big Bang Nucleosynthesis (O'Meara et al. 2001), and large- 
scale structure (Percival et al. 2001; Szalay et al. 2001), as 
well as the constraint on the size of mass fluctuations on the 
scale of the Jeans length at 2 < z < 4 imposed by the mea- 
sured flux per pixel distribution of the Lya forest (Meiksin 
et al. 2001). We use 512 3 particles on a 1024 3 mesh in a box 
of comoving side Lbox = 25/iT 1 Mpc. 

2.3 Sources of the UV background 

We consider two sources of the UV background: QSOs and 
stellar, such as from Lyman Break Galaxies (LBGs). 

We parametrize the QSO comoving luminosity function 
using the form introduced by Boyle et al. (1988), 

&(M,Z) = 10 o.4(/3 1 --l)[M*( z )-M] + 10 0.403 2 -l)[M*(z)-M] ' ( 5 ) 

where M*(z) represents an evolving break magnitude. The 
results of the 2dF QSO survey show that the QSO counts for 
z < 2.3 can be acceptably fit by a a pure luminosity evolu- 
tion model in which $* is independent of redshift, while the 
break magnitude evolves according to M*(z) = M*(0) — 
2.5(fci2 + k 2 z 2 ). For Q M = 0.3, Ov = 0.7, and h = 0.5, 
Boyle et al. (2000) obtain for the B— band luminosity func- 
tion /3i = 3.41, /3 2 = 1.58, $* = 0.36 x 10~ 6 Mpc~ 3 mag" 1 , 
Mb(0) = -22.65, fei = 1.36, and k 2 = -0.27. 

At 2 > 2.3, the QSO luminosity function is less well 
determined. We extend the luminosity function to higher 
redshifts using the results of the Sloan Digital Sky Survey 
for 2 > 3.5 (Fan et al. 2001a, 2001b, 2003). The survey 
limit is too bright to probe to fainter levels than the break, 
but shows a shallowing in the bright end of the luminosity 
function to /3i = 2.58, so that the pure luminosity evolu- 
tion model may no longer apply at the bright end. Since the 
pure luminosity model with f3\ = 3.4 is ruled out at only 
the 2-2.5 a level, and possibly only for the brightest QSOs 
(M1450 < -27) (Fan et al. 2001a, 2001b, 2003), we shall 
consider the pure luminosity evolution case as well. Since 
for the values of f3i and P2 we shall consider, eq. (5) pre- 
dicts that the UV background is dominated by sources with 
absolute magnitudes near or fainter than M* , the contri- 
bution of QSOs is poorly constrained by observations. We 
conservatively assume that the faint end of the luminosity 
function continues to follow pure luminosity evolution with 
P2 = 1.58, and fix the break luminosity M*(z) in eq. (5) to 
match the bright counts measured by the Sloan. 

The luminosity function reported by Fan et al. (2001a, 
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Figure 1. The comoving QSO luminosity function at z = 6 for 
the ACDM model with a bright end slope of @i = 2.58 (solid), 3.20 
(dashed), and 3.41 (dot-dashed), with the break magnitude M j 450 
adjusted to match the QSO counts at the bright end measured 
by the SDSS. 



2001b, 2002b) is in terms of the absolute magnitude M1450 
at 1450 A. To estimate the contribution of the QSOs to the 
UV background, we need to convert these values to lumi- 
nosities at the Lyman edge. We adopt the spectral shape of 
radio quiet QSOs as determined at lower redshift by mea- 
surements using the HST FOS (Zheng et al. 1997). For A < 
1050 A, we take /„ cx z/ -18 , while for 1050 A< A < 1450 A, 
we take /„ oc v~°" . The resulting break values Ml at the 
Lyman edge for /3i = 2.58 and 3.41 are given in Table 1. (The 
averages over the luminosity function have been performed 
over the range M1450 < —18.5, and the results adjusted to 
h = 0.7 for our ACDM model.) For (3i = 2.58, we find the 
proper emissivity decreases with redshift for z > 4 approx- 
imately as (1 + We also give values for the proper 
emissivity for /?i = 3.2. This is less than 2 a away from the 
fits of Fan et al. (2001a) and is able to match the emis- 
sivity required to reproduce the measured values of f a for 
4 < 2 < 6 using QSO sources alone. For this case, the proper 
emissivity evolves approximately as (1 + z) 1 ^ 4 for z > 4. We 
show the resulting comoving QSO luminosity functions at 
z = 6 for ft = 2.58, 3.20 and 3.41. 

We also consider a possible contribution from Lyman 
Break Galaxies or other stellar sources of Lyman contin- 
uum photons. Because UV continuum radiation shortward 
of the Lyman break is only observed from the brightest 
LBGs (Steidel et al. 2001), the full contribution of LBGs 
to the ionizing background is uncertain by a factor of a 
few for z ~ 3, and even more uncertain at higher redshifts. 
The spectral index of the emission, which affects the to- 
tal ionization rate, is also highly uncertain, depending on 
the processing of the continuum radiation through the in- 
terstellar medium of the galaxy and the IGM, and on the 



stellar populations of the galaxy. For the estimated ages 
of 10 8±0 B yr (Papovich, Dickinson & Ferguson 2001), the 
spectral shape is especially age-sensitive (Fioc & Rocca- 
Volmerange 1997; Bruzual 2001). For simplicity, we adopt 
a spectrally flat proper emissivity near the Lyman edge of 
eJ; BG k8x lO^ftergs^Hz^Mpc- 3 at z w 3 (Steidel 
et al. 2001), noting that this estimate was for an Einstein- 
deSitter universe, and so somewhat overestimates the contri- 
bution in the A-dominated universes we consider here. The 
uncertainty in the spectral index introduces an uncertainty 
in the photoionization rate of 30 — 50%. 



3 THE REQUIRED METAGALACTIC 
EMISSIVITY 

Within the context of any given set of cosmological param- 
eters, the flux distribution of the Lya forest predicted by 
numerical simulations is fully specified except for one un- 
certainty, the magnitude of the UV background responsible 
for the photoionization. Generally this is fixed by requiring 
the simulations to reproduce the measured average flux, or 
equivalently, the mean Lya optical depth f a . Usually this is 
done by specifying the required mean photoionization rate 
F per neutral hydrogen atom. In the presence of fluctuations 
in the UV background, however, the full distribution of F 
should be used, since f a depends non-linearly on F. This is 
particularly important when the distribution is broad. We 
predict the distribution below and will take it into account 
there. In this section, to simplify the calculation, we assume 
a constant value of F at any given redshift. 

The photoionization rate is related to the emissivity of 
the sources through 

where J„ is the angle- averaged intensity of the UV back- 
ground at frequency v, a v ~ <7l{v /vl)~ 3 is the photoion- 
ization cross-section, and vl is the frequency at the Lyman 
edge. The intensity J v is given by the proper emissivity e„ (z) 
of the sources according to (eg, Meiksin & Madau 1993), 

M*) = ^ £°° dz' (^±±) 3 e„ ( Z ')<exp [-r„(s, *')] )%W 

where (exp [—t v (z, z')]} is the average attenuation produced 
at v due to intervening photoelectric absorption by gas be- 
tween redshifts z and 2', dl p /dz' is the differential proper 
cosmological line element, v 1 = u(l + z')/(l + z), and it 
is assumed that the sources have turned on at some finite 
redshift z on . At high redshift in a flat universe, dl p /dz ~ 

(c/H )n-J /2 (i + zr 5 / 2 . 

Following Zuo (1992b) and Zuo & Phinney (1993), we 
introduce the attenuation length ro, defined by the ap- 
proximation (exp [— t v (z, z')]) w exp [— l p (z, z')/ro], where 
l p (z,z') is the proper distance between redshifts 2 and z' . 
Since our simulations do not include the effects of radiative 
transfer, the full contribution of Lyman Limit Systems is not 
accounted for. Their contribution to the total Lyman contin- 
uum optical depth from the IGM, however, becomes small 
compared with that from the Lya forest itself at the high 
redshifts we consider here (Meiksin & Madau 1993; Madau 
& Haardt 1996). As a test of their possible contribution, we 



© 2003 RAS, MNRAS 000, 000-000 



Constraints on the UV metagalactic emissivity using the Lya forest 5 



replaced the flux in pixels with tl > 1 by zero (fully black- 
ening them), and found this had a negligible effect on the 
values of ro obtained. 

It is informative to consider two limiting cases for the 
contribution of the sources to J„: 1. absorption limited (AL), 
for which the attenuation length is much smaller than the 
horizon, and 2. cosmologically limited (CL), for which the 
attenuation length much exceeds the horizon. For the AL 
case, the integral in eq. (7) reduces to J v = e v (z)r /4:ir. For 
the CL case, attenuation by the IGM is negligible, and the 
integral is limited by cosmological expansion and the age 
of the sources (avoiding Olber's paradox). In this case it 
is useful to represent the redshift evolution and frequency 
dependence of the emissivity as power laws: t v [z) = ex(l + 
z) v {i>/vL)~ as ■ Using the approximation above for dl p /dz, we 
then obtain J v (z) = e v (z)lv/4-w, where the effective proper 
distance in the universe over which sources contribute to 
J„(z) is 



*u = 



1/2 



(1+*) 



-3/ 



a !-[(! + *)/(!• 



i9/2+a s -77 



9/2 + as - T) 



.(8) 



Since the mean Lya optical depth constrains the photoion- 
ization rate V through the simulations rather than the emis- 
sivity directly, we re-express these results in terms of F. 
While the above approximations express the magnitude of 
the emissivity at the Lyman edge, they do not take into 
account the effect of attenuation at shorter wavelengths in 
case AL. The effect of intervening gas is to harden the spec- 
trum, since high energy photons are able to travel further 
than photons near the Lyman edge before they are photo- 
electrically absorbed (Haardt & Madau 1996). We account 
for the hardening by introducing the metagalactic spectral 
index omg, which replaces as for the AL case, but not for 
the CL case. The two limiting cases may then be combined 
into the approximation 

e L (z)a L 



r(z) 



(9) 



hp(3 + omg) 

where hp is the Planck constant, r e g = l/(l/ro +1//h), and 
'h = [(3 + «mg)/(3 + as)]lv- The cases AL and CL then 
correspond to the respective limits ro <JC Ih and ro ^> Zh- 

For any pair of ex and ro, T is specified through eq. (9), 
and therefore f a . Any pair that yields a desired f a , however, 
is not necessarily self-consistent. This is because any given 
value of T fixes not only f a but (exp(— tl)), and therefore 
ro, as well. It is thus necessary to solve self-consistently for 
values of ro as a function of ex- We show the dependence 
of ro and f a on ex for the results of our ACDM simulation 
in Fig 2. Here omg ~ has been taken to account for the 
hardening of the radiation field; the associated uncertainty 
in the amount of hardening, which will also depend on re- 
emission of Lyman continuum photons by the IGM itself, 
introduces an uncertainty of 30 — 50% in the estimates for 
ex - It is also assumed that the sources have turned on at very 
high redshifts, so that the limit z on — > oo may be taken in 
the expression for In. For small values of ex, the background 
intensity is attenuation limited (case AL). For large values, 
the intensity is limited by cosmological expansion (case CL). 
There is a rapid transition between the two regimes. For the 
constant comoving density source case, the transition ex- 
hibits critical behaviour, with a sudden transition from case 
AL to case CL as the emissivity is increased: a sufficiently 
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e L (10 28 h erg s" 1 Hz" 1 Mpc" 3 ) 




10 20 30 40 

e L (10 26 h erg s" 1 Hz" 1 Mpc" 3 ) 

Figure 2. The dependence of the attenuation length ro (top 
panel) and mean Lya optical depth f a (bottom panel) on the 
average emissivity at the Lyman edge for the ACDM model at 
z = 4, 5, and 6, as indicated. It is assumed that the (proper) emis- 
sivity evolves according to e„ = ex (1 + z) — 1 ("/"i) - 1-8 for z > 4, 
corresponding to QSOs (solid lines) or as e v = ex(l + z) A , cor- 
responding to a constant comoving source density (dashed lines). 
The horizontal bars in the upper panel mark the effective lumi- 
nosity horizons In of the sources (see text). 



large emissivity drives the attenuation length to such large 
values that attenuation becomes negligibly important. As 
attenuation becomes less significant, the mean Lya optical 
depth rapidly declines as well. Although the transition is 
less abrupt for the more gradual QSO evolution case, the 
ionization structure of the IGM is still very sensitive to the 
source emissivity. This sensitivity permits a tight bound to 
be set on the metagalactic emissivity. 

We show in Fig 3 the emissivity required to match the 
measured values of f a = 0.73, 2.1 and > 5.1 at the respective 
redshifts z ~ 4, 5 and 6. We use the results of Songaila & 
Cowie (2002) for the point at z = 4, estimating the error 
from the spread in their two fits to f a . The measured values 
at z ~ 5 and 6 are from Becker et al. 2001, with their quoted 
errors, and only the lcr lower limit is used at z ~ 6 based on 
Lya absorption. They set more stringent limits on f a at z ~ 
6 using Ly/3 and Ly7 absorption, but these are somewhat 
less firm due to absorption modelling uncertainties. 

The estimate of ex for Lyman Break Galaxies at z ~ 3 
(Steidel et al. 2001) exceeds the required value at z > 4 
by nearly an order of magnitude. If the proper emissivity 
were constant out to z = 4, attenuation by the IGM would 
become negligible, f a would be driven to very low values and 
the IGM would become nearly transparent to Lya photons. 
Either the proper emissivity of LBGs declines dramatically 
between z « 3 and 4, or the emissivity of LBGs at z « 3 has 
been overestimated. Ferguson, Dickinson & Papovich (2002) 
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Figure 3. (Top panel) The emissivity el at the Lyman edge 
required to match the measured mean Lya optical depth f a at 
z = 4, 5 and 6 (Becker et al. 2001; Songaila & Cowie 2002). The 
units of the emissivity are 10 26 /lergs" 1 Hz -1 Mpc -3 . The esti- 
mates based on the ACDM model are shown (points), with the 
error bars corresponding to the redshift ranges and measurement 
errors on f a . Only an upper limit is shown at z = 6, correspond- 
ing the the measured la lower limit on f a . The estimates for the 
emissivity from QSO sources alone are also shown, for /3i = 2.58 
(triangles), 3.20 (squares) and 3.41 (inverted triangles). (Bottom 
panel) The attenuation length ro, in proper Mpc (solid line) and 
the mean number No of QSO sources within the attenuation vol- 
ume (dashed line) predicted by the ACDM model are shown for 
f3i = 3.2, for which QSOs alone are adequate for providing the 
UV background at z > 4. 



argue that indeed the star formation rate in Lyman Break 
Galaxies declines sharply for z > 3. 

The estimated emissivities from QSOs assuming lumi- 
nosity slopes at the bright end of f3i = 2.58, 3.2 and 3.41 are 
also shown. While QSOs alone are not able to provide the 
required emissivity for (3i = 2.58, the preferred value from 
the Sloan (Fan et al. 2001a, 2001b), QSOs alone are easily 
adequate for f3i = 3.41, the value preferred by the 2dF at 
lower redshift. We find that the predicted QSO emissivity is 
very sensitive to (3i, scaling roughly as 6l ~ Pi' 5 - The QSO 
emissivity becomes sufficient for f}\ > 3.2 for z > 4 within 
the errors. The reason for the sensitivity to (3i is that as the 
luminosity function steepens, the low luminosity QSOs are 
able to provide an increasing proportion of the total emis- 
sivity. While the emissivity above Mi 450 varies little with 
/3i due to the constraints imposed by the number of QSOs 
detected at the bright end, the proportion of the total emis- 
sivity provided by QSOs fainter than Mi 460 increases from 
about 50% to nearly 80% as f}\ increases from 2.58 to 3.41. 
For (3i = 3.2, 75% of the total emissivity is provided by 
QSOs fainter than M* 450 . Thus most of the emissivity is 
derived from QSOs that are not directly detected. Yet it 
would seem unlikely that low luminosity QSOs, detected in 



lower redshift surveys, would disappear altogether. Our es- 
timates are based on the conservative assumption that the 
faint QSOs continue to follow pure luminosity evolution, al- 
though we have no direct evidence for this. 

The attenuation length ro is predicted to decay expo- 
nentially with redshift, r ~ 9 x 10 4 cxp [—1.5(1 + z)] Mpc 
for 4 < 2 < 6 and /3i = 3.2, as shown in the lower panel 
of Fig 3. The mean number No of QSOs within the at- 
tenuation volume similarly follows an exponential decay, 



iVo 



9 x 10 exp[— 4.35(1 + z)], with fewer than a 



gle QSO on average per attenuation volume by z = 6. If 
QSOs dominate the UV background, then large fluctuations 
in the background field and ionization level of the IGM will 
be introduced for z > 5. We next turn to the possible role 
played by a fluctuating background field on the absorption 
statistics of the IGM. 



4 UV BACKGROUND FLUCTUATIONS 

The stochastic nature of the discrete sources of the UV meta- 
galactic background will introduce local fluctuations in the 
intensity of the background. When these sources are suf- 
ficiently numerous, the fluctuations will be small. This is 
expected to be the case for any contribution from galaxies 
or other stellar sources (eg, intergalactic star clusters), be- 
cause of the large number density of these sources. Because 
of their rarity, however, the contribution from QSO sources 
will fluctuate at high redshifts. The degree of fluctuation is 
determined by the number of sources within a photoelectric 
attenuation volume, the region around the source to which 
Lyman continuum photons may propagate on average before 
being photoelectrically absorbed by the IGM. 

In order to compute the probability distribution /(J) of 
the UV background, we begin by computing the distribution 
/(J, N) of the UV background generated by an average of N 
sources within a finite sphere of radius R. Here J designates 
the angle averaged specific intensity of the radiation field 
at a given frequency and point in space, which we take to 
be the origin of our coordinate system. (We suppress the 
frequency subscript to simplify notation.) We work in the 
Euclidean limit. For any given realisation of N sources with 
specific luminosities Lk and distances rk, J is given by 



where j k = L k exp(-r fc )/(47rr fc ) 2 



(10) 

Tk/r is the op- 



tical depth due to photoelectric absorption by the IGM, 
and ro is the attenuation length. The distribution of J is 
given straightforwardly by the method of characteristic func- 
tions (Kendall & Stuart 1969) when the contribution of each 
source is statistically independent. The characteristic func- 
tion of f(J,N) is then (denoting an average over all the 
sources by (...)), 

f(t,N) = (cxp(iij)) (11) 
= J cxp(ifEk =1 j k ) [Hk =1 p(L k ,r k )dL k dr k ] ■ 

Here p(Lk,rk) dL k dr k denotes the probability of source k 
having a specific luminosity in the range (L fe , L k + dL k ) and 
lying within a sphere of radius R at a distance from the 
origin in the range (rk,rk + drk). It is given by 
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p(L k ,r k )dL k dr k 



4irri 



(4/3)^3 JW HL)dL 



dL k dr k , (12) 



and is identical for all the sources. The expression then sim- 
plifies to f(J,N) = I N (t,R), where, defining x — L/L* , 



I(t,R) = 



f B 3^2 rL nlax /L* 
JO ^ A min /L» 



dx(j>(x) exp(itj), 



(13) 



where we have used eq. (12) for p(L,r) (dropping the sub- 
script k), and 4>(x) is the normalized luminosity function of 
the sources 

3>(xL*)L* 



Six) — 



(14) 



where $(L) is the luminosity distribution per unit luminos- 
ity and unit (comoving) volume. 
Inverting f(J, N) we obtain 



fU. 



■■">-£/ 



exp(-itJ)I (t,R)dt. 



(15) 



Performing a Poisson average over N sources with mean 
N = (4n/3)R 3 n, where n = J ${L)dL, gives 



f(J,N) 



E 



N N e~ 
N\ 



-f(J,N) 



= ^- J e-x.p(-itJ)exp{N[I(t,R)-l]}dt.(16) 

Zuo (1992a) obtained eq. (16) using Markoff's method, 
although the derivation is more straightforward as shown 
here. He then presents solutions for the restricted case of an 
identical luminosity for all sources, no absorption within the 
attenuation volume (r k — 0), and R — ro- 

It is useful to recaste eq. (13) using the dimensionless 
variables s = tj* and r = r/ro, where J* = L*/(47rro) 2 . 
The expression for I(t,R) becomes 



I(s,R) 



3 iL d " 



/"^max ( 

= / dx6(x)l 



+ isx 



dx6(x)exp(isxT 2 e T ) 
dx<l>(x)^exp [isx(r /R) 2 e- R/r °] (17) 
-j=r J dr exp (is:ET~ 2 e~ T ) e~ T (2 + r) 



where x min = L min /L*, a; max = i max /L*, N = (A-k /?>)rln is 
the average number of sources within an attenuation volume, 
and we have performed an integration by parts to obtain the 
final expression. 

We may now take the limit R — + oo of eq. (17) to obtain 

lim N[I(s,R)-l] = isN / dxxS(x)G{sx) , (18) 

R— too J 

where we have defined the function 
CU)= I d^r 3 (0, 



f-OC 



(19) 



for which r is related implicitly to £ through £ = r~ 2 e~ T . 
Using eqs. (16) and (18), we obtain for the distribution func- 
tion of j — J I J* , 

f(j) = lim f(J,N)J* 

R. — 'oo 
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Figure 4. The probability density f(J) for the fluctuations in the 
intensity j for Nq = 1, 10 and 100 for the QSO luminosity function 
with fa = 2.58 and fa = 1-58 (solid lines). Also shown are the 
corresponding asymptotic limits f(j) ~ jNo{x 3 / 2 )j~ 5 / 2 (dashed 
lines). The distributions peak at lower for decreasing JVq. 



ds cos 



sN 



dxx6(x)Re G(sx) — sj 



x exp 



sNo 



dxxS(x)Im G(sx) 



(20) 



The leading asymptotic behaviours of the real and imag- 
inary parts of G(lu) for to <JC 1 are given by Re G(ui) ~ 

XT d ^ 3 (0 = 3 and Im G M ~ JT dtsm(^)C 3/2 = 
(2noj) 1/2 . These expressions may be used to derive the lead- 
ing asymptotic behaviour of f(j) for j 3> 1. We obtain 
f(J) ~ j N o(x 3/2 )j~ 5/2 (j > 1), found by rotating the inte- 
gral over s in eq. (20) into the complex plane and using the 
method of stationary phase. Here (x 3 ^ 2 ) is an average over 
the luminosity function 4>(x). The average value of j is found 
directly from eq. (16) to be (j) = dJ(J/J*)f(J,N) = 
3N () (x) (l -e- R/r °) -> 3iV (x) for R -> oo. Using the 
asymptotic form for f(j), we note that the rms fluctuations 
of j diverge at the bright end (j max > 1) like jUL- 

In this paper we solve eq. (20) numerically using the full 
QSO luminosity function as determined from observations. 
We show in Fig 4 the resulting distributions for No = 1, 10 
and 100, using eq. (5) with 0i = 2.58 and fa = 1.58. The 
peak of the distribution shifts toward lower intensity values 
as No decreases. 



5 THE Lya FOREST FLUX 
5.1 The pixel flux distribution 

The distribution function of flux per pixel has proven to 
be an extremely sensitive discriminator between Lya forest 
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Figure 5. Comparison of flux per pixel distributions (top panel) 
and cumulative flux distributions (bottom panel) at z = 4 (dot- 
dashed), 5 (dashed) and 6 (solid), normalized to the measured 
values of f a . The distributions sharply steepen at vanishing flux 
levels at high redshift. 



models. The number N of pixels available in a high reso- 
lution spectrum, like those from the Keck HIRES or VLT 
UVES, is on the order of 10 4 , so that the cumulative distri- 
bution of pixel fluxes may be measured to an rms accuracy 
of ~ N~ ' 6 « 0.01, although the finite width of the absorp- 
tion features tends to increase the rms somewhat because of 
the resulting statistical dependence between pixels (Meiksin 
et al. 2001). We show the pixel flux distributions from the 
simulations at z = 4, 5 and 6 in Fig 5. Because of the rapidly 
rising values of f a with redshift, the distributions become in- 
creasingly peaked at low flux values. 

Fluctuations in the UV background will affect the pixel 
flux distribution in two ways: a shift in its mean and a tilt 
in its shape. Matching the measured mean flux to the value 
predicted for a given set of cosmological parameters fixes 
the mean V for a simulation, so that neglecting fluctuations 
may lead to a mis-estimate of the mean T. 

At 2 = 4, the number of QSOs in an attenuation vol- 
ume is sufficiently large that the effects of background fluc- 
tuations will be negligible. We therefore only consider the 
higher redshift cases. We concentrate on the case f5\ = 3.2 
for the bright end of the QSO luminosity function, since the 
estimates in Table 1 suggests QSOs alone are able to meet 
the emissivity requirements to match the measured values 
of f a . 

We show in Fig 6 (upper panel) the required ionization 
rates, both with and without the effects of a fluctuating UV 
background. At z = 5, the fluctuations have only a small ef- 
fect on the required T, increasing it slightly above the case of 
smoothly distributed sources, for which T_i2 = 0.16. Once 
the effect of background fluctuations is included, it is no 
longer possible for QSOs alone to maintain the ionization of 



Figure 6. (Top panel) Estimates for the required mean 



H I photoionization rate r_i2, in units of 10 



required 



to match the measured f a at different redshifts. The lowest red- 
shift point is from Meiksin et al. (2001). The solid points are 
for an assumed homogeneous photoionization background. The 
open circle at z 5 takes into account fluctuations in the back- 
ground contributed by QSO sources for fi\ = 3.2, also including 
a small contribution from a smoothly distributed component. At 
z ~ 6, the lower upper limit is for the case of a homogeneous pho- 
toionization background. The limit marked Tg + Tqso includes 
both a smooth component and the effects of QSO fluctuations for 
Pi = 3.2. (Bottom panel) The measured f a at different redshifts. 
The lower points at z = 5 and 6 show the values of f a that would 
result assuming the same mean photoionization rate required by 
QSOs to match the measured values of f a , but neglecting the 
effect of the background fluctuations. The fluctuations substan- 
tially boost f a at z Ri 6. 



the IGM for /3i = 3.2; unless the metagalactic UV spectrum 
hardens to ckmg < 0, a small additional smooth component 
of Ts,-i2 = 0.025 must be added. The smooth component 
is an increase by 15% over the mean QSO ionization rate, 
a reasonable amount for the diffuse contribution expected 
to arise from the re-radiation of Lyman continuum pho- 
tons by the IGM itself. The difference is larger at z « 6. 
The total ionization rate must now increase from a mean 
of f _i2 = 0.057, to reproduce f a = 5.1, to 0.080, of which 
~ 40% is made up of a smoothly distributed component, 
again reasonable for re-radiation by the IGM. In the lower 
panel of Fig 6, a comparison is shown between the mea- 
sured values of the mean Lya optical depth and the values 
that would be obtained assuming the same required mean 
QSO rates as in the upper panel, but neglecting the role 
of the background fluctuations. The significant boost of the 
estimated mean Lya optical depth by the UV background 
fluctuations demonstrates the need for taking into account 
the role of the fluctuations at high redshifts (z £ 5). In par- 
ticular, this shows that care must be taken in ascribing an 
abrupt rise in the mean Lya optical depth to the reioniza- 
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Figure 7. Differences of the cumulative distributions of flux per 
pixel for the ACDM simulation with different assumptions for the 
origin of the photoionization background, at z = 5 (top panel) and 
z = 6 (bottom panel). The flux distributions are all normalized 
to the measured values of fa. The cumulative distribution differ- 
ences are with reference to the fiducial flux per pixel distribution 
for which a uniform photoionization background is assumed. The 
solid lines show the difference when fluctuations in a background 
dominated by QSOs are taken into account, adopting /3i = 3.2. 
The dashed lines show the difference assuming To = 10 4 K in- 
stead of To = 2 X 10 4 K for the temperature-density relation, and 
assuming a homogeneous UV background. 

tion epoch, as done by Becker et al. (2001), as such a rise 
may also be due in part to the UV background fluctuations 
resulting from a small number of sources in an attenuation 
volume. 

We note that the ionization rates in Fig 6 decay expo- 
nentially with redshift from z it 3.5, paralleling the decline 
in the number counts of QSO sources beyond this redshift. 

A comparison of the cumulative distributions of flux per 
pixel is shown in Fig 7. We show the differences in the cu- 
mulative distributions, as would be used in the Kolmogorov- 
Smirnov (KS) test. The fiducial models for z — 5 and z = 6 
assume a homogeneous UV background. We compare the 
homogeneous case with the case of /3i = 3.2, taking into ac- 
count both the fluctuations in the UV background from the 
QSOs and the small additional smooth component as found 
above. The effect on the cumulative distributions is small 
and beyond the current detection limits of the Keck and 
VLT, except possibly for large sample sizes. The differences, 
however, may be measurable with an ELT class telescope. 
We also compare the cumulative distributions (for smoothly 
distributed sources) for two assumptions for the IGM tem- 
perature. For the fiducial case, To = 2 x 10 4 K in eq. (4) was 
assumed. Decreasing To to 10 4 K significantly alters the cu- 
mulative distribution at z = 5, although has a smaller effect 
at z = 6. In principle, the uncertainty in the IGM temper- 
ature, which corresponds to a change in line-width of the 



Figure 8. The correlation function of the photoionization UV 
background generated by QSO sources alone assuming j3\ = 3.2, 
at z = 4 (dot-dashed), 5 (dashed) and 6 (solid). 

absorption features, may be corrected for by adjusting the 
line-widths to match the measured Doppler parameter dis- 
tribution or wavelet coefficient distribution (Meiksin 2000; 
Meiksin et al. 2001), so that this uncertainty may not be a 
limiting factor in detecting the effect of the fluctuations on 
the pixel flux distribution. We note that lowering the IGM 
temperature increases the radiative recombination rate and 
so increases the required ionization rate from T_i2 = 0.16 
to 0.22 at z = 5 and from r_i 2 = 0.057 to 0.069 at z = 6. 

5.2 Correlations in the UV photoionization 
background 

Zuo (1992b) showed that the discreteness of the sources of 
the UV background in the presence of attenuation will result 
in correlations between the background intensity at points 
separated in space. His result for the correlation function 
may be expressed as 

t ( \ 1 ( x2 ) r o Z* 00 , 1 , u + r/r _„ 

£ J v) = ^TTTT2 — / dw-log -j—e , (21) 

3N {x} 2 r J r/ro u u- r/r 

where the averages (. . .) are performed over the dimension- 
less QSO luminosity function 4>{x) (eq. [14]). In Fig 8, we 
show the correlation function for (3i = 3.2, for which nearly 
the entire UV background may be provided by QSOs. The 
correlation lengths at z — 4, 5 and 6 are 0.29, 2.2 and 
4.0 Mpc (proper), respectively, corresponding to the ve- 
locity differences (in the QSO restframe) of 125, 1250 and 
2850 kms -1 . The background correlations will induce corre- 
lations in the pixel fluxes as well, and so serve as a possible 
means for detecting the discreteness of the sources. The cor- 
relations will become somewhat diluted by any contribution 
from a more smoothly distributed component, as may arise 
from the diffuse background re-radiated by the IGM. We 
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treat the effect of UV background correlations in a separate 
paper (Meiksin & White 2003). 

5.3 Reionization of the IGM 

In this paper, we have considered QSOs as the possible 
sources that dominate the UV background at z Si 6. This 
does not require the QSOs to have initiated the ionization 
itself, however. The ratio of the Hubble time £h to the time 
for full radiative recombination in our ACDM model is 

te ™, 0.024 (^)(l + ^(^°' 65 , (22) 

where we have used eq. (4) and approximated the 
radiative recombination coefficient as cka ~ 4 x 
10- 13 (T/10 4 K)-°- 7 cm 3 s _1 . At z = 6, the ratio exceeds 
unity only for p/p Si 3.5. At the mean IGM density, the ra- 
tio is less than unity for z < 11. If an early generation of 
AGN or stars pre-ionized the IGM, for instance, it would be 
possible for the generation of QSOs at z < 6 to maintain the 
ionization before the IGM had time to fully recombine. It is 
none the less interesting to address the question of whether 
these QSOs would have been sufficient to have done the 
reionization as well. The number density of ionizing pho- 
tons generated by the QSOs over a Hubble time, assuming 
/3i = 3.2 and t\ = 1/4, isn 7 («) = 6.8x 10" 4 (l+z)" 5/4 cm~ 3 . 
Comparing this with the total number density of hydrogen 
atoms «h(z) gives 

lh^^ 4000 (^Y\ l + z) -iy\ ( 2 3) 

n H (z) ^ 0.020 y v ' v ' 

This crosses unity at z = 6.04. If these QSOs were respon- 
sible for the reionization of the IGM, then they only just 
could have done so at the highest redshifts to which they 
have been seen. In the presence of clumping, reionization by 
the QSOs becomes problematic. If, however, the IGM is still 
partially ionized at z ~ 6 following an earlier epoch of reion- 
ization, then the QSOs would be adequate for renewing the 
ionization of the IGM, as the required number of ionizing 
photons scales like the neutral fraction of hydrogen atoms. 



6 SUMMARY 

We have investigated the dependence of the ionization struc- 
ture of the IGM at high redshifts (z Si 4) on the emissivity 
of the sources of the UV photoionization background in the 
context of a ACDM universe. We show that the ionization 
structure is highly sensitive to the emissivity of the sources 
at the Lyman continuum edge. For low emissivity values, the 
intensity of the UV background is limited by the attenuation 
of Lyman continuum photons by the IGM. For large values 
of the emissivity, the IGM becomes essentially transparent 
to Lyman continuum photons, and the UV background is 
limited instead by cosmological expansion and the age of the 
sources. We show there is a rapid transition between these 
two extremes with increasing emissivity. If the UV back- 
ground is dominated by Lyman Break Galaxies at z Si 4, 
then their proper emissivity must dramatically decline from 
the value estimated at z w 3 (Steidel et al. 2001), or they 
would render the IGM nearly transparent at the Lyman edge 



and result in a much smaller mean Lya optical depth than 
measured. 

We have also explored the possibility that the UV back- 
ground is dominated by QSO sources at high redshift. We 
show that within the observational errors on the high red- 
shift QSO luminosity function, QSOs may provide a sub- 
stantial component to the UV background, and may even 
dominate it. If they do dominate, then their sparseness will 
result in substantial fluctuations in the UV background. We 
derive the distribution function for the fluctuating inten- 
sity in an infinite universe including the effects of attenua- 
tion, and show that the distribution significantly broadens 
for z Si 5. The fluctuations in the UV background will in- 
crease the estimated Lya optical depth over the value that 
would be estimated assuming a homogeneous background. 
This increases the demand placed on the ionization rate re- 
quired to reproduce a measured mean Lya optical depth 
and should be taken into account when comparing with es- 
timates of the QSO emissivity. It also shows that care must 
be taken in attributing a sudden rise in the Lya optical 
depth to the reionization epoch (Becker et al. 2001), since 
a sudden rise may result instead from a rapidly decreasing 
attenuation length and the resulting UV background fluctu- 
ations associated with the reduced number of sources that 
dominate the background. There is a secondary effect of the 
fluctuations on the pixel flux distribution of the Lya forest, 
resulting in a small distortion. Although the effect is too 
small to be measured in a single QSO spectrum from ex- 
isting facilities like the Keck or VLT, it may be detectable 
by combining many spectra, or may need to wait until the 
arrival of a future generation of extremely large telescopes. 

The large mean Lya optical depth measured by the 
Sloan at z ~ 6 (Becker et al. 2001) requires fewer than a sin- 
gle QSO within an attenuation volume on average. This is a 
significant measurement in that it implies that if QSOs dom- 
inate the UV background, there will be large fluctuations 
in the mean Lya optical depth from region to region, over 
spatial scales of ~ 4 Mpc, the expected correlation length 
of the background radiation field produced by the QSOs at 
2Ri6. This offers a possible means of distinguishing between 
models for the origin of the UV background dominated by 
either rare sources like QSOs, or by relatively more copious 
sources like galaxies. We pursue this topic in a companion 
paper (Meiksin & White 2003). 

We have also shown that if QSOs dominate the UV 
background at z « 6, they would just have had sufficient 
time to provide an adequate number of Lyman continuum 
photons to ionize the hydrogen, but not at larger redshifts. 
The estimate, however, is based on the minimal counting 
requirement of the production of one ionizing photon per 
hydrogen atom, and neglects the possible role of clumping 
of the gas. On the other hand, simulations show that much 
of the volume of the IGM will be underdense and so more 
readily photoionized. If the IGM were ionized at an earlier 
phase and had not yet fully recombined by z ~ 6, then even 
in the presence of clumping, QSOs may be adequate for re- 
newing the ionization of the IGM. We also note that the 
UV background fluctuations we expect, if QSOs dominate 
the ionization of hydrogen, will apply to helium ionization 
as well, for which it is even more certain that QSOs were the 
sources of reionization. A complete description of the reion- 
ization of the IGM will require accounting for the nature 
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of the ionizing sources, their immediate environments, and 
the inhomogeneous structure of the IGM. Reviews of recent 
work in this evolving subject are provided by Madau (2000) 
and Loeb & Barkana (2001). 

In order to match the statistics of observed spectra in 
detail, it will ultimately be necessary to incorporate radia- 
tive transfer into the simulations both to allow for the ex- 
pected spatial correlations in the UV ionizing background as 
well as to reproduce the correct gas temperatures (and line 
widths) , particularly if helium is completely photoionized at 
only moderate redshifts. Several other physical effects may 
also be important, such as heating sources in addition to 
QSOs, or galactic feedback like winds. 



z 




M L 










Pi = 


2.58 


3.20 


3.41 


2.58 


3.20 


3.41 


4 


-22.86 


-23.98 


-24.23 


3.23 


7.96 


9.77 


5 


-22.11 


-23.44 


-23.73 


2.68 


8.19 


10.54 


6 


-21.52 


-23.02 


-23.35 


2.36 


8.62 


11.52 



Table 1. Break luminosity M£ at the Lyman edge and proper 
emissivity at the Lyman edge required to match the Sloan QSO 
counts at different redshifts, assuming a QSO luminosity func- 
tion slope at the bright end of f3\ = 2.58, 3.20 and 3.41 and 
a homogeneous UV background. The units of the emissivity are 
lO^hergs-iHz^Mpc- 3 . 
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